Dynamic spin-Hall effect and driven spin helix for linear spin-orbit interactions 
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We derive boundary conditions for the electrically induced spin accumulation in a finite, disordered 
2D semiconductor channel. While for DC electric fields these boundary conditions select spatially 
constant spin profiles equivalent to a vanishing spin-Hall effect, we show that an in-plane ac electric 
field results in a non-zero ac spin-Hall effect, i.e., it generates a spatially non-uniform out-of-plane 
polarization even for linear intrinsic spin-orbit interactions. Analyzing different geometries in [001] 
and [110]-grown quantum wells, we find that although this out-of-plane polarization is typically 
confined to within a few spin-orbit lengths from the channel edges, it is also possible to generate 
spatially oscillating spin profiles which extend over the whole channel. The latter is due to the 
excitation of a driven spin-helix mode in the transverse direction of the channel. We show that while 
finite frequencies suppress this mode, it can be amplified by a magnetic field tuned to resonance 
with the frequency of the electric field. In this case, finite size effects at equal strengths of Rashba- 
and Dresselhaus SOI lead to an enhancement of the magnitude of this helix mode. We comment on 
the relation between spin currents and boundary conditions. 

PACS numbers: 72.25.Dc 85.75.-d, 75.80.-|-q 



I. INTRODUCTION 

Electron systems with spin-orbit interaction show a va- 
riety of spin-electric effects arising from the coupling be- 
tween (orbital) charge and spin degrees of freedom. The 
most prominent examples are the spin-Hall effectii^i^ 
and current induced spin polarization^ii, both of which 
have received substantial interest due to their potential 
to generate and control spin polarization with electric 
fields. This type of electrical control is a prerequisite for 
integrating spin effects into standard lithographic semi- 
conductor structures and, ultimately, utilizing the spin 
degree of freedom as a carrier of information^ 

The spin-Hall effect (SHE) manifests itself 
experimentallj*2i^i^ as current induced spin polar- 
ization (CISP) at the edges of a Hall-bar (in the 
absence of a magnetic field). Initial theoretical studies 
of the SHBiS in 2D electron systems have focused on 
linear intrinsic Rashba- and/or Dresselhaus spin-orbit 
interaction (SOI) and interpreted this boundary spin 
accumulation in terms of a spin current^ (defined as 
a symmetrized product of spin and current densities) 
flowing transverse to the applied electric field. However, 
these arguments have been plagued by ambiguities, such 
as equilibrium spin currentsii and the absence of spin 
conservatio n^^i^^ in systems with intrinsic SOI. Explicit 
diagrammatic calculationsi^ii^ii^ for disordered systems 
and a more general, non-perturbative argumcn t^^i^'^i^^ 
show that the spin-current is absent in systems with 
standard linear-in-momentum SOL— 

A more straightforward approach is to calculate the 
quantity directly measured in experiments: the spatially 
and time resolved spin densityJ^i^ii^ In weakly disor- 
dered systems with Ep t~^,Aso (where Ep is the 
Fermi energy, r the momentum relaxation time, and 
Aso the spin-orbit splitting) the spin density is described 



by spin diffusion equations derived in Keldysh2£i2ii^ 
or density matrix approaches J^i22ii^ These equations 
have been used to study various effects, such as the re- 
sponse to an electromagnetic wave^^ spin currentSf2£ 
spin relaxation^i^ boundary spin accumulation for 
(Jq17j^j19j28^ and abruptly switched^Si^i electric fields, 
and more general interface problemsi^i^i 

A significant difference between charge and spin diffu- 
sion, as described by these equations, is the existence of 
spatially oscillating spin density modes. For instance, a 
gradient of the out-of-plane spin density acts as a torque 
on the in-plane spin and vice versa, leading to a peri- 
odic spatial modulation of both in- and out-of-plane spin 
densities with a period given by the spin-orbit length 
Aso- General solutions of the spin diffusion equations 
are damped spatial spin density oscillations with a pe- 
riod given by the spin-orbit relaxation length Ago. An 
example of such periodic modes in diffusive systems was 
first described in Ref.ll^ (see, in particular, Eq. (7) there) 
for the case of equal strengths of the Rashba and linear 
Dresselhaus SOI. For this particular case and in the ab- 
sence of the cubic SOI, these modes are long-lived and 
static and are thus referred to as persistent spin helixi^ 
Modes of this type have recently been observed 

However, when analyzing these equations for a specific 
geometry, e.g., in a narrow channel for the case of the 
SHE, the weight of these oscillatory modes in the solu- 
tion is determined by boundary conditions (BCs). For in- 
stance, assuming vanishing polarization at the boundary 
one obtains an oscillatory behavior of the spin density;^ 
resembling the spin profile measured in Ref. dj. On 
the other hand, for a von-Neumann boundary condi- 
tion (vanishing normal gradient of the polarization), the 
spin profile is spatially uniform. Thus, the existence of 
the SHE depends crucially on the BCs. This circum- 
stance motivated a number of studies where BCs for sys- 



terns with SOI were derived microscopically, both in the 
difFusiveiiii^i^i^^i^ and ballistio^ regimes. 

It has been show n^^i^^ that BCs (for hard-wall spin- 
conserving boundaries) in disordered^ systems with lin- 
ear SOI and for dc electric fields require the spin den- 
sity to be equal to its value in the bulk, i.e., far away 
from the boundary, and, thus, lead to a spatially uni- 
form spin profile. This null result is consistent with zero 
spin currents ji^ii^ii^ii^i^S The experimentally observed 
dc spin accumulation^ in 2DEGs thus requires an expla- 
nation accounting for both extrinsio^i^ and (cubio^^) 
intrinsic effects. That a spin current is finite at finite 
frequencies and for linear SOIsr^i^i^ however, hints at 
the presence of boundary spin accumulation in ac solu- 
tions. In this article, we focus on the intrinsic mechanism, 
and show that a dynamic SHE, i.e., boundary spin po- 
larization induced by an ac voltage, is present even in a 
minimal intrinsic model. 

The dynamic SHE arises due to the excitation of spa- 
tially non-uniform spin diffusion modes. In the Hall-bar 
geometry, these modes are excited by a spatially uniform 
ac electric field and lead to accumulation and spatial 
oscillations of the spin density close to the boundaries. 
Analyzing these modes as a function of SOI strengths 
and in the presence of an external, in-plane magnetic 
field, we find a spin diffusion mode which is a finite- 
frequency analog of the persistent spin helix i^ii^i^i^^ 
The relaxation length of this mode -while finite for 
generic linear SOIs- becomes infinite when the Rashba 
and Dresselhaus SOI strengths are equal and when the 
magnetic field is tuned to resonance with the frequency 
of the electric field. This particularly robust mode, 
originating from electric-dipole-induced spin resonance 
(EDSR)^>ia^'^>iiii^>i^>ii>i^ gives rise to a spatially os- 
cillating spin profile which extends infinitely far away 
from the Hall-bar boundary. This driven spin helix has 
the same spatial oscillation period as the persistent spin 
heli3*2ii22iS but, whereas the latter is static, the former 
oscillates in time at the frequency of the applied bias. 
The prediction of a driven spin helix is one of the main 
results of this paper. 

Using a linear response approach, we solve the prob- 
lem of a hard-wall boundary in a disordered 2D electron 
gas in the presence of an ac electric field. The deriva- 
tion of the BCs is similar to the one in Refs. [l^ and 
[29I . We find that while the bulk polarization is reduced 
at finite frequencies, the BCs require the polarization at 
the boundary to have a larger value. The spin polariza- 
tion is, thus, no longer spatially uniform: there is a spin 
accumulation at the boundary and spatial oscillations de- 
caying towards the bulk of the sample. The amplitude 
of this spatial oscillations at frequency lo is proportional 
to w/r, where T is the Dyakonov-Perel spin relaxation 
rate. Since typically T <C t~^, where r is the transport 
time, the dynamic SHE becomes pronounced even for 
frequencies ojt <C 1. Analyzing different geometries and 
SOIs, we find that it is possible to excite a predominantly 
oscillatory mode for equal strength of the Rashba and 
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FIG. 1: Left: Conducting channelinfinite in the e^ij -direction 
and of width L in the 6^;^ direction, ac Electric field E(ti;)||ea;2 
induces boundary spin accumulation. An external magnetic 
field bo, applied parallel to E, gives rise to EDSR (see 
Sec. IVII[) . Right, (a): a "standard" [001]-grown quantum 
well with the [110] crystal axis taken along the a; 2 -direct ion. 
The bulk polarization f2(eEr) oc ea;^(a + /3) points along 6^;^ 
[cf. Eqs. (dl and (O]. Right, (b): a [110]-grown quantum 
well with E||[ilO] along e^^- The internal field [see Eq. Q] 
fl(eEr) has both in-plane (due to the Rashba SOI) and out- 
of-plane (due to the Dresselhaus SOI) components. 



Dresselhaus SOI -a driven spin helix described above. 

The paper is organized as follows. In Sec. |TT1 we in- 
troduce our model and formulate the linear response for- 
malism for SHE. In Sec. IIIIl wc sketch the derivation of 
the integral equation for the spin density, which is then 
used to derive the diffusion equation and its boundary 
conditions. (A more detailed derivation is deferred to 
appendix [XI) In Sec. llVi we derive boundary conditions 
in the presence of ac electric field and comment on the 
relation between spin currents and these boundary con- 
ditions in Sec. |Vl In Sec. IVII we calculate the spatially 
resolved spin profiles at finite frequencies in various ge- 
ometries in [DDI]- and [110]-grown quantum wells. Gener- 
ation of a driven spin helix under the conditions of EDSR 
is discussed in Sec. IVIII . 



II. PRELIMINARIES 

We consider a disordered 2DEG confined to a quantum 
well (QW) channel of width L (sec Fig. [1]) with non- 
interacting electrons of mass m and charge e. The system 
is described by the Hamiltonian 



Zm 



bo 



V. 



(1) 



Here, p = {pi,p2,0) is the in-plane momentum, r2(p)i = 
^ijPj is a linear, vector-valued function of p describing 
spin-orbit interaction, 2bo = g^B{Bi, B2,0) is a mag- 
netic field (equal in magnitude to the Zeeman energy) 
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applied parallel to the 2DEG, and <t = {a^,(j'^,a^) are 
the Pauli matrices (and ct° = 1). The disorder potential 
V due to static short-ranged impurities randomly dis- 
tributed over the channel is characterized by the mean 
free path I = rpp/m, where r is the scattering time and 
Pi? is the Fermi momentum. 

We calculate the impurity-averaged, spatially- 
dependent spin density 5"(r) = cr*(5(r — x) due to 
in-plane ac electric field E(a;) — Eo(5(ti-' — ojo)- As it 
will be shown below, the overall magnitude of S is 
determined by the bulk spin polarization due to CISP 
far away from the boundary. We therefore briefly discuss 
CISP in different geometries. We define the nominal 
polarization 

Sfc EE -jy2rj(eE(cj)r) , (2) 

(with J/ = m/2TT being the density of states per spin) 
which at zero frequency {ujq ~ 0) coincides with the bulk 
polarization!^ In this case, Sf, is simply a paramagnetic 
spin response to an effective magnetic field O(eEoT). The 
latter is the internal field due to the electrically induced 
drift momentum eEr and SOI. 

Both the magnitude and direction of Sb depend on the 
SOI mechanism. We consider two cases (see Fig.[T|): the 
"standard" [001]- and [110]-grown QW. The Rashba SOI 
(with strength a) due to an asymmetry in the confine- 
ment potential has the same form in both cases and is 
assumed to be tunable. The Dresselhaus induced fields 
J^£>,[ooi], i^D,[iio] are in-plane and out-of-plane in the 
[001] and [110]-grown QWs, respectively. For conve- 
nience, we define = 2apFT, S^p = 2Ppft as the ratios 
of the mean free path and spin precession length due the 
Rashba and Dresselhaus SOIs, respectively. The vector 
couplings of the SOIs are described by 

/ a + (3 0\ 
^^[001] - -(«-/?) (3) 
\ 0/ 

for case (a) in Fig. [1] and 

/O a 0\ 

^[110] = -a (4) 
\0 13 J 

for case (b). In case (a), the bulk polarization Sb oc 
—e^i {a + P) points along the (negative) xi-axis. When 
the Rashba- and Dresselhaus SOIs are of comparable 
strength, i.e., a ~ +f3 (or a ~ — /?), constructive (de- 
structive) interference between the two SOI mechanisms 
occurs.— In this case, one spin component [along xi (X2)] 
becomes conserved. A similar situation occurs in the 
[110]-grown QW, where the out-of-plane spin is conserved 
if the Rashba SOI is relatively small. Here the bulk polar- 
ization points out-of-plane and is, thus, easier accessible 
in optical measurements . i^s 

The induced spin density S'^ (r) is described by coupled 
spin diffusion cquation a^^i^"'^'^ which can be derived in 
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FIG. 2: (color online) Characteristic wave numbers Q of the 
homogeneous solutions S5(ri) = ss.oe*^"^^ and sa = S6(ri) = 
S6,oe*'^''^ of Eq. ([7| as a function of a//3 for ldt = 10"^ and 
= 0.1 in a [001]-grown quantum well. For a = —f) (in- 
dicated by arrow), the wave numbers have small real parts 
|Re6'5| ~ |Re\/— 2iLjr| <^ 1, which implies a nearly undamped 
oscillatory mode. Under EDSR conditions, Re^s vanishes 
identically at a-|-/3 = [cf. Eq. (|32|) ]. The presence of modes 
with almost imaginary wave numbers leads to an oscillating 
spin profile, as shown in Fig. [3] 

the Keldysh22iiSii2^ or density matrix formalisms F^SiS^ii^ 
As a starting point for the derivation of the bound- 
ary conditions, we present here an alternative derivation 
based on a diagrammatic linear response approach. The 
detained derivation is deferred to Appendix El We ob- 
tain an integral equation for the spin density 

S\v)- Sl^iujTSl + j <fxX'^{v,Ti){S^{-K)- si) (5) 

valid in the regime l/EpT <^ 1, where 

^^^(r,x) = ^tr{a'^Gt+Jr,x)a'^Gt(x,r)} . (6) 

Here, Ep is the Fermi energy, tr{. . . } denotes the trace 
over spin s, Vj = ^ + ilkjO''' is the velocity operator con- 
taining a spin-dependent term, and G^^^ the impurity- 
averaged, retarded/advanced Green functions at energy 
E. Note that for w = the integral equation ([5]) depends 
only on the combination — S;, so that the spatially uni- 
form solution S' = Sf, is immediate. The uniform spin 
profile is equivalent to the absence of the SHE, whose 
presence would cause a spatial modulation of S at the 
boundary. 

III. DIFFUSION EQUATION 

Far from the sample boundary, the impurity-averaged 
Green functions and, hence, the kernel X'"'(r, x') « 
g-|r-x \/i decay on the scale of the mean free 

path I, which is the shortest length scale of the diffu- 
sion problem. The behavior of S on scales larger than I 
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can, thus, be found by expanding: S(x) w S{r) + (x — 
r)idriS{r) + ^{ii — r)k{x — r)idr^dr,S{r). In this way, one 
obtains the coupled spin diffusion equation 

[-tLu + T-DAr]iS{r)-Sb) (7) 
- 2[b - ppnilVr)] X (Sir) - Sb) = iujSb , 

where D = vpl/2 is the diffusion constant and F*-' = 
[tr{(rif2^)}(5,;j — {ilH'^)ij]2p'^T is the spin relaxation ten- 
sor. 

We now apply Eq. ([7]) to the two specific geometries 
in Fig. [T] (a) and (b). Assuming translational invari- 
ancc along ex2, we find for the [001]-grown QW with 
E||[110]||e., 

[-IUJ - Ddl^ + F_] - Sb) + C-dr,S'^ = lujSb , (8) 
[-IUJ - Ddl^ + F+] 52 = , (9) 
[-iuj - DSif.^ + F+ + F_] S"^ - C-dr,S^ = , (10) 

where F± 2p%T{a ± (3f, C± = 2pfI{P ± a), ujl = 2bo 
and Sb = -2veET{a + /3). 

In case (b) of a [110]-grown QW with bo, E| | [110] | |e^, 
we find 
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FIG. 3: (color online) Real (solid line) and imaginary parts 
(dashed line) of the out-of-plane spin density S'\x)/Si so- 
lution of Eqs. (O-dlOl), CSI] in the standard QW (Fig. [T] 
a) are shown for lot = 10~*, = 0.1, ^/j = —0.095, and 
L — 500L As Rashba- and Dresselhaus SOI interfere destruc- 
tively for this case, the bulk polarization is much smaller than 
the value discussed in Sec. IVII all other parameters being 
equal, Sb = 0.005 fj,m~^ instead of 1.1 fj,m~^ . Inset: In-plane 
polarization S^{x) (along the internal field i7(eEr)) in the 
same situation. 



[^Dd^,^-zu; + r[ + r',]is' -Si) 
-[C'.dr, + u;l + V^2]is' - s!) 
[-Dd^^-iuj + T[ + r'^]s^ = 0, 

[-Ddl~iu; + 2T',]{S'~S',) 

+[c^dr, +LUL- Vr^Ks^ - sl) 



iuSl 



b ' 



iujSf 



b ' 



(11) 

(12) 
(13) 



where Sb = -2veET{a,Q, (3), F'^ = 2p|,r/32, F!, = 
2p|,ra^, C[ = PfW, C2 = Ppla, and ujl = 26o- Note 
that in this geometry the Dresselhaus SOI adds to Sb, 
whereas for E||[001] the electric field does not couple to 
the Dresselhaus term4^ 



IV. BOUNDARY CONDITIONS 

The diffusion equation ([7|) has to be supplemented with 
boundary conditions. These match the bulk solutions of 
the diffusion equation ([7]) with the solution of the in- 
tegral equation Eq. ([5]) in the region 1/pf <C xi ^ / 
close to the boundary. Here, we follow the approach 
used in Refs. [l^ and [2^. We choose a;i = as the 
boundary and construct the impurity-averaged Green 
functions G'^/"*(x,x') which satisfy the Dyson equation 
(x|[(u; -Ho- Hso - S)G]|x') = 5{x - x') with Hq be- 
ing the Hamiltonian in the absence of SOI and S be- 
ing the self-energy due to impurity scattering42i^ We, 
moreover, impose the hard-wall, spin-conserving bound- 
ary conditions C?(x, x')|^j ,j.'^=o = for either argument 
at the boundary. 

To 0th order in the SOI, these conditions are satisfied 



by image constructions Gq^^ 



G 



R/A 
b,0 



G'^J^, where 



bulk and Gj^^^'^(x, x') = G^^^{yi, {—x'i^x'2)) is the Green 
function mirror-reflected at the boundary. Neglecting 
Friedel oscillations of the self-energy at the boundary, 
which fall off as 1/ ^ppxi, the Green functions G^^^ 
constructed in this way satisfy the Dyson equation to 
leading order in 1/EpT. 

To 1st order in Hso the Green functions is found 
as Gi = GqHsoGo- By construction, G(x,x') = 
[Go + Gi](x, x') satisfies the boundary conditions and the 
Schroedinger equation to linear order in the spin-orbit in- 
teraction. Performing a Fourier transform of the Green 
bmction G(x,x') = /(ip2G(.Ti,x'Jp2)e*P^("^^-^2)/(27r) 
along the boundary, we find 



Gq^^{xi,x[\p2) = 



(14) 



where p^ = yj2'm{E ± i/2t — p\/2m) with p2 being the 
momentum along the channel. To first order in Hqq, we 
find 



G^^^{xi,x[\p2) 



(X] 



~x',) 



(15) 



Gjj Q is the impurity-averaged Green function in the 



where the dots stand for additional terms that do not 
contribute to the integrals below since they are odd in 
the longitudinal momentum p2. 

We are now in a position to derive boundary conditions 
using the Green functions from Eqs. (fTil) and We 
take the limit r ^ of Eq. Q and expand S"(x) k, 
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S*(r) + {xj - rj)^5'(r) in the integrand. This yields 

= iLurSb + iB- 1)(5(0) - Sb) + C,—S{0), (16) 

orj 



where the coefficients 



dx[dx'2X'"'{x,:ic')\a 



(17) 



da;;da;'2X'"'(x,x')(x' - a;j)U^o (18) 



are obtained from the spin-spin correlation function X^"^ 
in Eq. © evaluated with the Green's functions satisfying 
the boundary conditions. In symbolic notations, X cx 

G'b,o^b,o + ^hfi^bfi ~ ^b,o^b,o ~ ^bfi^hfi- that 
for Lo = the diffusion equation Eq. ^ and the boundary 
conditions Eq. (jl6p depend only on the combination 5' — 
Sfc, so that the spatially constant solution S = Sf, is 
immediate. In particular, there is no spin accumulation 
close to the boundary in that case in agreement with the 
literature on the linear intrinsic spin-Hall effect. ^'^i^^i^^i^^ 
When calculating the spin-spin correlation function 
X^'^one encounters mixed terms of the form GG* , which 
oscillate as a function of xi with a period of 1/pf- To 
determine S on length scales larger than I, we neglect 
these oscillations. This way, we find the BCs 

WaS' = -2pFTn{n)mermj [s' - , (19) 

where n is a unit vector normal to the boundary and 
where we have neglected terms proportional to ujt <C 1. 



V. SPIN CURRENT 

In this section, we show that a definition of the spin 
current in terms of an SU(2)-covariant derivative is con- 
sistent with both the boundary conditions and the dif- 
fusion equation. This definition is equivalent to defining 
the spin current as the commutator - in contrast to the 
conventionally used anficommutator - of spin and ve- 
locity. To see this, we define a Hermitian spin current 
operator as follows 

J^{S) = -I?,^''(r) = -i[mvso,u ^"(r)] (20) 

where we have introduced the covariant derivative^ 
= d/dxi-_ — i[A, J with the non-abelian gauge po- 
tential Ai = —milkia'' and S^{r) = (T^6{r — x) is the 
spin density operator. [Note that Jl' differs by a factor 
of mass m from the conventionally defined product of 
velocity and spin.] From Eq. ()20p we obtain a spin cur- 
rent J by replacing iS by S{r) — Sb in the second fine of 
Eq. jSni), i.e.. 



Jl'(r) = -I?r' (^'''(r)-^;'' 



(21) 



where -2?,'''' = S'T^'d/dn + 2ml]feiefc,,r,'- The BCs in 
Eq. are then equivalent to the requirement that the 
normal component of J vanishes at the boundary, i.e., 
n-J"(r)|n-o = 0,77=1,2, 3. 

Using the definition, Eq. ([2T|) . one finds that both the 
diffusion equation [for this see also RefUBl, Eq. (I7|), and 
the boundary conditions, Eq. (jl9p . can be written in 
terms of the covariant derivative as 



^5" + DVT J = 



n - J" 



= 0. 



(22) 
(23) 



Thus, spin diffusion with linear SOI has a (formal) anal- 
ogy to charge diffusion: In charge diffusion, both the 
diffusion equation p = DVj for the charge density p and 
the BCs n • j = contain the same charge current j. The 
current j = Vp is given in terms of the spatial derivative 
of the density. Analogously, the spin current is given as 
the SU(2)-covariant derivative of S^. 

In Ref. [29l . in an attempt to identify a spin current di- 
rectly from the diffusion equation, Eq. ([7| was rewritten 
(for bo = 0) in the form 

-iLuS'^ + [T{S - Sb)]" - DV • J" = , (24) 
where the "spin current" 

Jl' = -^S^r) - AmnMekn'^S"^' (r), (25) 

however, differs from by a relative factor of 
2. This discrepancy is resolved when the definitions 
Eqs. ((20)) . ((2T|) are used making the introduction of two 
different spin currents J and J unnecessary. 



VI. SOLUTIONS OF THE DIFFUSION 
EQUATION 

First; we obtain a solution of Eq. ^ in an infi- 
nite sample. In this case the bulk Green's functions 
g5^(x,x') 



G^^^(x — x') are translationally invariant 
and; thus, A^'' = J (Px'X'^'^{r,x') becomes independent 
of r. The spatially uniform ansatz 



5oc = [r - ' 



(26) 



solves both the integral equation [Eq. ^ and the diffu- 
sion equation [Eq. The same result for the polar- 
ization at finite oj was found in Ref. using a kinetic 
equation and in Refs. [stIIssI in the linear response for- 
malism. Remarkably, Soo is not simply given by the ac 
internal field n{pd) corresponding to ac drift momen- 
tum pd = eE(ijj)T/(l — iuor), but depends on the spin 
relaxation rate. Therefore, the deviation of Soo from 
Sh becomes appreciable already at a relatively small fre- 
quency ~ F^ip rather than at a much higher frequency 
Lo ~ r^^, which marks the dispersion of p^. Note also 
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FIG. 4; Real(upper panel) and imaginary parts (lower panel) 
of the out-of- plane spin polarization S'^{x)/Sb for frequencies 
LOT = 0.1,0.3,0.5,0.7,0.9, 1.1 xlO"'^ (grey to black lines), = 
0.1, = 0.003, L = 2800Z, and Si, = l/^m"^. The geometry 
with E||[ll0] is shown in Fig. [T](b). Horizontal dashed lines 
mark ac bulk polarization according to Eq. (|26p for the same 
parameters. 



that there is no dc bulk polarization at \a\ — \f3\, if the 
limit of a — !■ ±/3 is taken before the limit of lu ^ [see 
Ref. [U for a discussion of this point] . 

We now estimate the magnitude of Sb. We choose pa- 
rameters similar to the (low-mobility) sample employed 
in Ref. |2 except for a higher mobility and a lower sheet 
density. With a = 1.0 x 10~^^ eVm, sheet density 
712 = 1.0 X 10^^ m~^, and transport mean free path r = 
5 X 10~^^ s and choosing E = 5 mV//xm, we obtain the 
bulk polarization due to Rashba SOI Sb.a = iveEra = 
l.l^m^^, or about 1 spin per ^m^ ( Sb/n2 = 0.1%). The 
magnitude of Soo and, as we will see below, the magni- 
tude of the spatially non-uniform terms in the solution 
are proportional to Sb- Depending on the geometry and 
on whether the Rashba and Dresselhaus SOIs add con- 
structively or destructively, the overall amplitude of the 
spin oscillations and edge spin accumulation is modified. 
In case (a) in Fig. [U one finds Sb = Sb^a[^ + /3[ooi]/'^) 
while in case (b) Sb = Sb,a(X,^, ot), where /9[ooi] 

and /3[iio] is the Dresselhaus SOI strength in the [001]- 
and [110]-grown QW, respectively. 

We now focus on the position-dependent spin profile in 
a semiconductor channel of finite width. As before, we 
assume translational invariance along the channel so that 
the diffusion equation Eq. ([7]) becomes an inhomogeneous 
ordinary differential equation 



in the transverse coordinate ri , where the differential op- 
erator L{dri) is defined by Eqs. ([7]) and ([?7|) . The solu- 
tion 



-c*=Sfc(r) 



(28) 



consists of the uniform part S'oo, given by Eq. 
(inhomogeneous solution), and a linear combination of 
k = 1,2,..., 6 cigenmodcs = Sk^e^''^'^, satisfying 



X(Vr)s(r) = 0. The wave numbers 
order) in case (a) are given by^ 



.,6 



H,2 



^3,4,(5,6) 



+ (-)! 



2D 



ioj) 

2D (r+ - 2iuj) - c'i 



(in arbitrary 



(29) 



+ (-) ip 2\ -2D (r+ - 2w) CI 



D-^V\ 



L(S,J[5(ri)-Sfc] ^JwrSfa 



(27) 



Some of 6'fe are shown in Fig. [2] as functions of a//?. The 
real and imaginary part of the wave number arc respon- 
sible for exponentially growing (decaying) and oscilla- 
tory parts of the mode, respectively. The coefficients 
are determined by the boundary conditions in the 
form Mc = - {{B - \){S^ - Sb),-{B - 1){S^ - Sb)) 
where M is a 6 x 6 -matrix obtained by inserting the 
general solution into Eq. (fTH]) [see also Eq. (jB3P in Ap- 
pendix |B] . The coefficients c'^ determine the magnitude 
of the non-uniform part of S, i.e., if all c'^ are zero the so- 
lution is spatially uniform. Although explicit expressions 
for c are too lengthy to be displayed here, the scaling of 
c with iu can be found on general grounds. Indeed, all 
the entries of the matrix Af~^diag((i3 — 1), — (B — 1)) are 
of order 1. The order of magnitude of Ck is thus given 
by \Soo — Sb\ ^ {u)/r)\Sb\, where the latter holds for 
w ^ r. The non-uniform part of S (proportional to the 
c's), thus, scales linearly with oj for uj -^T and becomes 
appreciable at the frequency scale w ~ F ^ r~^. 

A solution for S in a [001]-grown QW (Fig. [1] a) is 
shown in Fig. O The electric field E is along the [110] 
axis and the strengths of the Rashba and Dresselhaus 
SOIs are chosen as a — /3, so that the wave numbers (cf. 
Fig. ^ are almost imaginary. In this case, oscillations of 
the out-of-plane spin density S^ extend almost over the 
entire channel. Simultaneously with a approaching —f3, 
however, the internal field ri(eEr) cx a + /3 and, thus, the 
overall amplitude Sb of the spin density becomes small. 
In other words, suppression of the damping rate Re^,; oc 
\a + (3\ close to the special point a = —(3 competes with 
a suppression of the overall amplitude, so that a purely 
oscillatory mode cannot be excited in this geometry. 

Figure|4]depicts the polarization profile in a wide [110]- 
grown QW as shown in Fig.[T](b), where the bulk polar- 
ization (due to the Dresselhaus SOI) is out-of-plane. For 
a weak Rashba SOI, the wave numbers of the character- 
istic modes are almost real, i.e., the modes are strongly 
damped. As a result, the polarization close to the bound- 
ary is substantially larger than the bulk value given by 
Eq. 
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FIG. 5: Polarization in the EDSR geometry E, B]|y for case 
shown in Fig. [2a) with Sb — 0.02/^m~^. Upper panel; 
Im5''^(a:: = 201) (black) and ImS^{x = 201) (grey curve) are 
shown as a function of ujlt. Resonance is seen at ztcoLT ~ 
LOT = 10"^. Parameters of the [001]-grown QW: ia = -0.08, 
= 0.1, L = 100/. Lower Panel: Density plot of Re5'^(2;), 
lmS^{x) as a function of x and for the same param- 

eters. 



VII. EDSR AND DRIVEN SPIN HELIX 



We now focus on electric-dipolc- induced spin resonance 
(EDSR}^iia>ii'i^>^^'^>^ii^ in the finite Hall bar geom- 
etry. We calculate the spin polarization S due to a simul- 
taneous effect of ac electric field and dc magnetic field bo, 
both along the channel. The directions of the fields are 
chosen in such a way so that the internal field r2(eE(a;)T) 
and bo are perpendicular This geometry is suitable for 
an observation of electrically driven Rabi oscillations of 
the spin polarization between the directions along and 
opposite to bo-^i^sdo 

We focus on case (a) in Fig. [1] The magnetic field 
bo leads to an equilibrium polarization (Pauli paramag- 
netism) Sbo cc X2<j-'L = 2bo in the longitudinal direction 
of the channel. In addition, the polarization in the bulk 
of the sample (transverse to bo) is modified. In the ge- 
ometry of Fig. [1] (a) with bollea;^, we find for the bulk 



r_(r+ + r_ 



Sb 



-iuj)) 



(30) 



-uj^ -iLu{r+ + 2r^) + T^_ 



r,r_ 



where r± = 2p|,r(a ± /3)^. In the absence of the mag- 
netic field, i.e., for ujl = 0, Eq. ([50]) reduces to Eq. ((^ . 
Additionally, the characteristic modes change due to bo. 
The wave numbers 9 are determined by the requirement 
of vanishing eigenvalues 



1 



1 



4{ujL-0C^)' 



= 
(31) 



of the diflFerential operator L{9) defined by Eqs. ([7]) . ([27ll . 

We focus on the case of a w — /3. Expanding to first 
order in r+/ {ujl~ 9C-) <C 1, one finds 



^^1,2 =T/"V2r(r+ 



iC- T V2V£' [r+ -2i{uj + cjl) 
2D 



75,6 



2D 



(32) 



At resonance, i.e., for uj]^ = lj, the wave numbers 9^^ ^ 
in Eq. (|32p become purely imaginary because r+ — 
for a — —j3. The modes 55^5 E^re thus completely un- 
damped oscillations of the spin density with wave length 
= l/2m(/3 - a) [cf. Fig [2]. Note that in the consid- 
ered case of a = —[3 the Hamiltonian commutes with the 
longitudinal spin ([i/, cr^] = 0), i.e., the [/(l)-symmetry 
described in Ref. l32l remains intact; however, the SU{2)- 
symmetry used in Ref. [U to demonstrate the existence 
of the persistent spin helix is broken in the presence of 
bo. 

Figure [5] shows a profile of the spin polarization under 
EDSR conditions. At resonance (w^ = iw), the overall 
amplitude of the out-of-plane polarization is enhanced. 
This enhancement becomes particularly strong for r_|- « 
occurring at a = — /3. 

Solving the difi^usion equation [Eq. ([7])] to first order 
in r+ for the case a w — /3, u « -l-w^, we obtain the 
following expression for the spin density close to reso- 
nance [uj « +i-ul) 



Sin) 



{Sl-Sl-iSi,){^p~^a,) 
sinh{LR/l)R 

-^ri/Xso cosh{R{L-ri) /I) 



g»(L-ri)/A30 C0sh(i?ri//) 



(33) 
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where R = ■\/rr+ — 2i{uj — lol)t. Equation ([33]) de- 
scribes a spin density wave along the transverse direction 
of the Hall bar with wave length X^q and an amplitude 
proportional to 1/ smh{LR/l) R. We discuss this result 
in more detail below. Inserting Eq. ((30)) for 5'oo into 
Eq. ([33]) , setting oj = lul, and expanding the hyperbolic 
functions in Eq. (j33|) for a narrow channel with width 
L <|; XgQ — l/2m{a + (3), one obtains the dominant 
a-dependence of 5' around the a w — /3 point as 



a + (3 



(34) 



wherein = [— icj/r_](— 2z^e£'r)(f;3— ^q,)Z/L depends only 
on the combination (3 — a. Here, we introduced a phe- 
nomenological linewidth r'°" = r;^'=" + ^T^^"" + 0{{ujTf) 
to model the regularization of the amplitude of at 
a-|-/3 = 0, which for V"'' = would diverge as l/(a + /3). 
For a + f3 = 0, the relaxation mechanisms due to linear 
intrinsic SOIs, which are dominant for generic a ^ ±/3, 
are ineffective, and finite spin relaxation rates and 
r^'^'' of the x\ and xi spin components, respectively, are 
due to an extrinsic or cubic Dresselhaus SOI. 

Equation describes a spin density wave S^ir-C) at 
frequency oj with a spatial profile of the form e~^'^^l^so . 
The real and imaginary parts of have stationary nodes 
separated by the shortest of the two SO lengths, i.e., \~^q- 
In addition, the spin profile is subject to a quantization 
condition: S is proportional to a factor 1 — e'^^^so, which 
vanishes for L = iirNX^Q (with N being an integer) and 
becomes maximal for L ~ (2iV+ l)7rAgQ. The profile de- 
scribed by Eq. ([M)) arises due to an excitation of the spin 
helix modes 55^6 under the EDSR conditions. The spa- 
tial oscillations of these modes have the same "magic" 
wave number 9 = as the static persistent spin 

helix j^^i^^ However, whereas the persistent spin helix is 
time- independent, the spin profiles in Eq. (P51) . (|34p oscil- 
late also in time at each point ri with the frequency of 
the applied electric field. [The explicit time-dependence, 
e.g., S{ri,t) cx sin {ti/X^q + cuot) for L = {2N + l)TTXgQ 
and for E(t) = Eg cos{LUot), is obtained by inverse Fourier 
transform of Eqs. (|33p . ((34)l .]. This driven spin helix is a 
generalization of a static spin helix structure to the time- 
dependent case. 

Spatial quantization due to the Hall-bar boundaries, 
moreover, leads to further enhancement of the amplitude 
of the spin helix modes in the EDSR regime. The ampli- 
tude (a + P)/ [i^a + C/3)^ + 2rr'<='^] is infinite for a^-P 
in a model with strictly linear SOI, i.e. for T'^'^^ = 0, but is 
regularized by the next-to-leading order effects due to cu- 
bic Dresselhaus and extrinsic SOIs, giving rise to a finite 
linewidth p''^''.— Such an enhancement of the amplitude 
of the driven spin helix close to the a = — /3 point in rel- 
atively narrow QWs may be observable, e.g., by optical 
techniques i^i^ 



VIII. CONCLUSIONS 

In conclusion, we have described several signatures 
of electrically induced spin polarization and the spin- 
Hall effect due to linear spin-orbit interactions. We 
have shown that the spin-Hall effect and edge spin 
accumulation- while being absent for dc electric fields- 
becomes finite for time-dependent electric fields. In par- 
ticular, we have found that boundary effects can extend 
over the whole sample due to driven spin helix modes for 
the case of the linear Rashba and Dresselhaus spin-orbit 
interaction being of equal strengths. The amplitude of 
these helix modes as a function of the spin-orbit inter- 
action strengths is strongly enhanced due to due spa- 
tial quantization under the conditions of electric-dipole- 
induced spin resonance. 
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APPENDIX A: SPIN DIFFUSION EQUATION 

We start from the impurity averaged Kubo formula 
(for EpT ^ 1) for the spin density 





ZTT 



— / d^x' 



5"'6{r - x') 
d2a;2mTA'''(r,x)i:'^''(x,x')]7"'(x'), (Al) 

where solid lines denote impurity averaged Green's 
functions G^^^ and dashed lines denote correla- 
tors of impurity potentiali^'^'^ The first term 
of Eq. ([XT]) is the "bubble" diagram -y^ir) = 
tr{(r|(T''G|^^_^%G';|^|r)} Ej{lj), i.e., a spin response to 
the electric field in the absence of vertex corrections. The 
latter are decribed the diffuson _D^''(r,x), which is de- 
fined by the integral equation 



L'^''(r,x') ^ 



5^'''5{Y - x') 
2mT 



d?vX^P{v,y)DP^{y,^), 
(A2) 



where X^^^ is given by Eq. ([6]). Iterating Eq. (|Aip once 
with the help of Eq. ([A2|) . we find 



S\v) 



e 

2^ 



(2TOr) j d^x'D'''{r,ic')-f''{ic'). (A3) 
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Multiplying Eq. (jA2p by ^2mr7'^(x') and integrating 
over x', we obtain the integral equation for the spin den- 
sity 



S\r) = — y(r) + J d2;^x-(r,x)5'^(x) , (A4) 

which can be further simplified by partially evaluating 
the "bubble " term 7* (r) in Eq. (|A4p . We define the spin- 
momentum correlation functions^ 

f |!^tr(a"G«(r,x)-j^G^(x,r)l (A5) 



2niT 

1 — iiUT 



mdxi 



-(1 + iUT)fl 



V3 ' 



(A6) 



where Eq. (jA6[) is obtained by evaluating Eq. jASj) us- 
ing the bulk Green's functions of an infinite sample. In- 
serting the definition of the velocity operator Vj = ^ + 
^kja'', we obtain (e/27r)7'(r) = - J d'^xX'''{r,x)Sl; + 
27nT{e/2'K)Y'^r)Ej. We can rewrite Eq. (|^ as 



+ [Y{r) - Y^f 2veEjT . (A7) 



From now on, we treat the regions close to the boundary 
and in the bulk separately. In the bulk, one obviously 
has \Y{v) — Yb] = and arrives thus at Eq. ([5]). At the 



boundary, the Green's functions Gq 



R/A 



G 



R/A 
b,0 



G 



-R/A 
b,0 



constructed in Sec. |V] have to be used to evaluate 7*(r), 
y^^(r). Neglecting terms oscillating with a period of 
1/pF, as described in Sec. IIV| one finds that y*-'(r)£'j = 
Y^-' Ej to linear order in the SOI. Therefore, the last term 



in Eq. (jA7P vanishes. Consequently, Eq. \K1\ turns into 
Eq. ^ and can be used for the derivation of both the 
bulk diffusion equation and the boundary conditions. 



APPENDIX B: BOUNDARY CONDITIONS 

For the coefficients B and G in Eq. describing a 
boundary with normal vector n, we found 



m^^'^in) = [B- ir = —2pFTn^{n)e^,^, (Bl) 

TT 

Gf"(n) = (B2) 



where we neglected terms proportional to lot <^ 1. We 
define a 6 x 6 matrix 



M 



(5B(n) + eiC(n))si,oe''i' 



(5S(-n) + eiC(-n))si.oe''i 



(5B(n) +e2C(n))s2,oe''2'' 
I 



{SB(A) + eeCin))se,oe''e'' 



(5B(-n) + e6C(-n))s6,oe''e' 



(B3) 



and a vector A = (Ao,Al), where Aq^l 
6B{±h){S CO — S&). Inserting the general solution S 



Soo + CfcSo,fce '^''^ into Eq. p6)) . the BCs can be rewritten 
as Mc = -A. 



1 M. I. D'yakonov and M. I. Perel, JETP Lett. 13, 467 
(1971). 

2 Y. K. Kato, R. C. Myers, A. C. Gossard, and D. D. 
Awschalom, Science 306, 1910 (2004). 

^ V. Sih, R. C. Myers, Y. K. Kato, W. H. Lau, A. C. Gossard, 
and D. D. Awschalom, Nature Physics 1, 31 (2005). 

* H.-A. Engel, E. I. Rashba, and B. I. Halperin, in Handbook 
of Magnetism and Advanced Magnetic Materials (John Wi- 
ley & Sons Ltd, 2007). 

^ L. S. Levitov, Y. V. Nazarov, and G. M. Eliashberg, Zh. 
Eksp. Teor. Fiz. 88, 229 (1985). 

® V. M. Edelstein, Solid State Comm. 73, 233 (1990). 
Y. K. Kato, R. C. Myers, A. C. Gossard, and D. D. 
Awschalom, Appl. Phys. Lett. 87, 022503 (2005). 

® D. D. Awschalom, D. Loss, and N. Samarth, eds., Semicon- 
ductor Spintronics and Quantum Computation (Springer, 
Berlin, 2002). 

^ J. Wunderlich, B. Kaestner, J. Sinova, and T. Jungwirth, 



Phys. Rev. Lett. 94, 047204 (2005). 
^" J. Sinova, D. Culcer, Q. Niu, N. A. Sinitsyn, T. Jungwirth, 

and A. H. MacDonald, Phys. Rev. Lett. 92, 126603 (2004). 
" E. I. Rashba, Phys. Rev. B 68, 241315(R) (2003). 

S. 1. Erlingsson, J. Schliemann, and D. Loss, Phys. Rev. B 

71, 035319 (2005). 

0. Chalaev and D. Loss, Phys. Rev. B 71, 245318 (2005). 
J. Inoue, G. E. W. Bauer, and L. W. Molenkamp, Phys. 
Rev. B 70, 041303(R) (2004). 

R. Raimondi and P. Schwab, Phys. Rev. B 71, 033311 
(2005). 

1^ O. V. Dimitrova, Phys. Rev. B 71, 245327 (2005). 

1. Adagideli and G. E. W. Bauer, Phys. Rev. Lett. 95, 
256602 (2005). 

A. G. Mal'shukov, L. Y. Wang, C. S. Chu, and K. A. Chao, 
Phys. Rev. Lett. 95, 146601 (2005). 

V. M. Galitski, A. A. Burkov, and S. DasSarma, Phys. 
Rev. B 74, 115331 (2006). 



10 



^° E. G. Mishchenko, A. V. Shytov, and B. I. Halperin, Phys. 
Rev. Lett. 93, 226602 (2004). 

B. A. Bernevig, J. Orenstein, and S.-C. Zhang, Phys. Rev. 
Lett. 97, 236601 (2006). 

Y. Tserkovnyak, B. I. Halperin, A. A. Kovalev, and 
A. Brataas, Phys. Rev. B 76, 085319 (2007). 
A. A. Burkov, A. S. Nunez, and A. H. MacDonald, Phys. 
Rev. B 70, 155308 (2004). 

I. Adagideh, M. Scheid, M. Wimmer, G. E. W. Bauer, and 

K. Richter, N. J. Phys. 9, 382 (2007). 

A. Shnirman and I. Martin, Europhys. Lett. 78, 27001 

(2007). 

Y. V. Pershin, Physica E 27, 77 (2005). 

T. D. Stanescu and V. Galitski, Phys. Rev. B 75, 125307 

(2007). 

E. I. Rashba, Physica E 34, 31 (2006). 

O. Bleibaum, Phys. Rev. B 74, 113309 (2006). 

R. Raimondi, C. Gorini, M. Dzierzawa, and P. Schwab, 

Sohd State Comm. 144, 524 (2007). 

N. P. Stern, D. W. Steuermann, S. Mack, A. C. Gossard, 
and D. D. Awschalom, Nature Physics 4, 843 (2008). 
J. Schhemann, J. C. Egues, and D. Loss, Phys. Rev. Lett. 
90, 146801 (2003). 

J. D. Koralek, C. P. Weber, J. Orenstein, B. A. Bernevig, 
S.-C. Zhang, S. Mack, and D. D. Awschalom, Nature 458, 
610 (2009). 

V. A. Zyuzin, P. G. Silvestrov, and E. G. Mishchenko, 
Phys. Rev. Lett. 99, 106601 (2007). 

H.-A. Engel, B. I. Halperin, and E. I. Rashba, Phys. Rev. 
Lett. 95, 166605 (2005). 
^® R. Raimondi and P. Schwab, Eur. Phys. Lett. 87, 37008 
(2009). 

M. Duckheim and D. Loss, Nature Physics 2, 195 (2006). 
M. Duckheim and D. Loss, Phys. Rev. B 75, 201305(R) 
(2007). 

R. L. BeU, Phys. Rev. Lett. 9, 52 (1962). 
*° E. I. Rashba and A. L. Efros, App. Phys. Lett. 83, 5295 
(2003). 

''^ E. I. Rashba and V. L Sheka, in Landau level spectroscopy, 
edited by G. Landwehr and E. L Rashba (North-Holland, 
Amsterdam, 1991), p. 131. 

*2 Y. K. Kato, R. C. Myers, A. C. Gossard, and D. D. 
Awschalom, Nature 427, 50 (2004). 

"^^ M. Schuhe, J. G. S. Lok, G. Denninger, and W. Dietsche, 



Phys. Rev. Lett. 94, 137601 (2005). 

Z. Wilamowski, H. Malissa, F. SchafHer, and W. Jantsch, 
Phys. Rev. Lett. 98, 187203 (2007). 

L. Meier, G. Sails, L Shorubalko, E. Gini, S. Schon, and 

K. Ensslin, Nature Phys. 3, 650 (2007). 

R. Winkler, Spin-orbit Coupling Effects in Two- 

Dimensional Electron and Hole Systems (Springer, 2003). 

E. Akkermans and G. Montambaux, Mesoscopic physics of 

electrons and photons (Cambridge University Press, 2007). 

J. Rammer and H. Smith, Rev. Mod. Phys. 58, 323 (1986). 
'^^ I. V. Tokatly, Phys. Rev. Lett. 101, 106601 (2008). 

O. E. Raichev, Phys. Rev. B 75, 205340 (2007). 
""^ M. Duckheim, D. Loss, M. Scheid, K. Richter, Inang 

Adagideli, and P. Jacquod, (unpublished) (2009). 

0. Chalaev and D. Loss, Phys. Rev. B 80, 035305 (2009). 
^'■^ M. Lee, M. O. Hachiya, E. Bernardes, J. C. Egues, and 

D. Loss, Phys. Rev. B 80, 155314 (2009). 
'^^ M. Duckheim and D. Loss, Phys. Rev. Lett. 101, 226602 
(2008). 

J. H. Bardarson, I. Adagideli, and P. Jacquod, Phys. Rev. 
Lett. 98, 196601 (2007). 

1. Adagideli, J. H. Bardarson, and P. Jacquod, J. Phys. 21, 
155503 (2009). 

P. G. Silvestrov, V. A. Zyuzin, and E. G. Mishchenko, 
Phys. Rev. Lett. 102, 196802 (2009). 
''^ In quantum wells with more than one subband the spin- 
Hall current can be nonzero. See ReflH. 
Qualitatively different spin effects occur in systems with 
orbital phase coherence, i.e., in mesoscopic disordered^ or 
ballistic systems^i^ and on length scales smaller than the 
mean free path at ballistic boundaries.— 
Below, we focus on the behavior of S for a « +/3 {a ~ —f)). 
In a [110]-grown quantum well, similar behavior can be 
found when the Rashba SOI strength is close to zero and 
the electric field is taken along the [110]-direction (the 
[001] -direction). Therefore, we discuss case (a) in more de- 
tail. 

We also studied the Hanle geometry when n(eE(Lj)r)j|bo; 
however, the effect of bo becomes appreciable only for fields 
as large as 6o ~ in that case. 
®^ A similar regularization from effects outside the linear 
model occurs for the lifetime of the static spin heli5*2i. See 
Ref . [23. 



